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Abstract. - We show that the glass transition predicted by the Mode-Coupling Theory (mct) 
is a critical phenomenon with a diverging length and time scale associated to the cooperativity 
of the dynamics. We obtain the scaling exponents i/ and z that relate space and time scales 
to the distance from criticality, as well as the scaling form of the critical four-point correlation 
function. However, both these predictions and other well known MCT results are mean-field in 
nature and are thus expected to change below the upper critical dimension dc — 6, as suggested 
by different forms of the Ginzburg criterion. 



One of the most striking property of glassy materials is the extremely fast rise of their 
relaxation time (or viscosity) as the temperature is lowered or the density increased [1] . The 
basic mechanism for this spectacular slowing down is not well understood, but it is reasonable 
to think that it is intimately related to cooperative effects. The dynamics becomes sluggish 
because larger and larger regions of the material have to move simultaneously to allow for a 
substantial motion of individual particles. Although this qualitative idea has pervaded the 
glass literature for many years [1] , it is only quite recently that a proper measure of cooperativ- 
ity (and of the size of the rearranging regions) was proposed and measured experimentally [2] 
and in numerical simulations [3-5]. The idea is to measure how the 'unlocking' events are cor- 
related in space; technically, this involves a four-point density correlation function (see below) 
from which one can extract a growing dynamical correlation length [4-6]. Recent extensive 
numerical simulations in Lennard-Jones systems have confirmed the crucial importance of this 
growing length scale for the dynamics of the system [4-6]. Furthermore, the four-point cor- 
relation function is found to have scaling properties similar to those expected near a critical 
point [4,5,7], suggesting that the physics of the glass transition should be understood as a 
critical phenomenon dominated by large scale fluctuations. This ingredient in fact appears 
in various forms in several recent phenomenological approaches [7-10]. This is, at first sight, 
in plain contradiction with the Mode-Coupling Theory (mct) of the glass transition, mct 
is considered by many to be the closest to a first principle, microscopic theory of glasses yet 

© EDP Sciences 



2 



EUROPHYSICS LETTERS 



achieved, with many quahtative and quantitative successes in explaining various experimen- 
tal and numerical results [11-13]. However, freezing in mct was argued repeatedly by its 
founders to be a small scale phenomenon, the self-consistent blocking of the particles in their 
local cages. Since no small wavevector singularities seem to exist in mct, power-laws and 
scaling are only expected in time but not in space [12]. This is surprising since on general 
grounds a diverging relax;ation time can only arise from processes involving an infinite number 
of particles. This point of view was challenged by Pranz and Parisi [14] in the context of the so- 
called 'schematic' mct (sec also [8] for an early insight). This simplified version of the theory 
is (formally) equivalent to describing mean-field spin glasses with three body interactions, for 
which the physics is known in great details. At the (crgodicity breaking) critical temperature 
Tc, the curvature of the relevant TAP states is known to vanish [15]. Therefore, one expects, 
and indeed finds, that a susceptibility diverges when T ^ Tc- This susceptibility turns out 
to be the precise analogue, in the context of this spin model, of the four-point correlation 
function mentioned above. By analogy with usual second order phase transitions, the analysis 
of [8, 14] suggests that the mct freezing transition is in fact accompanied by the divergence 
of the correlation length of the four-point correlation function. 

In this letter we show that the MCT dynamical transition in finite dimensions must indeed 
be understood as a critical phenomenon: dynamical correlations become long-range both in 
time and space, in agreement with the insight of [8,14]. We obtain the MCT dynamical scaling 
exponents v and z that relate space and time scales to \T — Tc\, as well as the scaling form 
of the critical four-point correlation function. However, these results (as well as all other 
quantitative mct predictions) are mean field in nature and change in dimensions less than 
dc = 6 due to long-wavelength fluctuations, as confirmed by different forms of the Ginzburg 
criterion. Our strategy is similar to the one used for ordinary critical phenomena. Consider 
the ferromagnetic Ising transition as an example. In that case one can show that there exists 
a certain functional of the magnetisation field such that (a) its first derivative leads to exact 
equations for the magnetization and (b) its second derivative is the inverse of the spin-spin 
correlation function. In general one cannot compute this functional exactly but one can 
guess its form using symmetry arguments, or compute it approximately in a diagrammatic 
expansion. Its simplest version corresponds to the Ginzburg-Landau free energy functional. 
The saddle point equations for the magnetization then leads to the mean field description of 
the transition. One finds in particular that the singular behaviour at the transition is related 
to the vanishing of a 'mass'. This has two important implications: (1) the spin susceptibility 
diverges at the transition, (2) because of the vanishing mass the corrections to the mean 
field result computed by adding more diagrams blow up whenever d < dc = ^, i.e. spatial 
fluctuations change the critical behavior. We shall show that exactly the same scenario takes 
place within the mct of the glass transition, except that the order parameter is now a two point 
function, the dynamical density- density correlation function. The analogue of the spin-spin 
correlation function is therefore the following four point correlation: 

Gi{r,t;lT) = {p{0,0)p{S,T)p{f,t)p{f+5,t + T)) - (p(0,0)p(5,r))(p(f,t)p(f+5,i + r)), (1) 

where p{x, t) is the density fluctuation at position x and time t and (•) denotes an average over 
the dynamics. Its intuitive interpretation (for example in the case where (5 = 0) is as follows: if 
at point an event has occurred that leads to a decorrelation of the local density over the time 
scale r, what is the probability that a similar event has occurred a distance r away, within the 
same time interval r, but shifted by t7 In other words, 6*4 (r, t: S, r) measures the cooperativity 
of the dynamics. Different values of r and t allow one to measure the full space-time structure 
of dynamical cooperativity, and are the analogue of the space and time separation entering the 
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Fig. 1 - Basic vertex of the MCT (a), and a ladder diagram that contribute to the divergence of the 
four-point susceptibihty (c). Other ladders can be obtained by combining the eight elementary blocks 
shown in (b). The full lines represent correlation functions, whereas the mixed dashed-fuU lines are 
response functions. 

standard spin-spin correlation function in critical dynamics. The quantities S and r, on the 
other hand, serve to define the 'order parameter', i.e. the density-density correlation function 
C{5,t) — (p(0, 0)p((5, r)). Previous studies have focused on the case t = 0, and S smaller 
than the particle radius, to which one can associate a 'susceptibility' Xii^) by integrating 
over space. More generally, we shall focus on a wave-vector dependent, a.c. susceptibility as: 

The quantities G4 or X4 can be computed as in the ferromagnetic case inverting the second 
derivative of an appropriate functional. Within any field theoretical derivation of MOT (e.g. 
the Das and Mazenko formulation [16] or more heuristic derivations [13,17]), the functional 
alluded to above can be constructed as [18]: 

F{G) - -^TrlogG + ^TtG-'G + <S>2pi{G) (3) 

where G and Go are compact notations for the full and the bare propagator of the theory, 
and <&2P/(G) is the sum of all two particle irreducible Feynman diagrams (that cannot be 
decomposed in two disjoint pieces by cutting two lines) constructed with the vertices of the 
theory and using the full propagator as line. The first derivative of F gives back the exact 
equations on G, whereas the four point function is obtained inverting its second derivative, 
which is nothing else than G~^G~^ (coming from the derivative of the first term in (jJJ) minus 
the derivative of the self-energy with respect to the two point function. Thus, since MOT is 
tantamount to only retaining the 'bubble' diagram, the four point functions are obtained from 
ladder diagrams, such as the one drawn in Fig. WC')- 

In order to understand the mechanism that leads to a divergence of the four-point correla- 
tion, consider the ladder diagram shown in Fig. 1. The n*^ order contribution to X4(fc, oj, K, il) 
reads: 

A.Uk,^Xn) = ld%d%.^d^k^,ld.^^^^^^^ 

■^fe.i^(^l' fcri-i;t^l,t^ri-l)A^fc.^(fcn-l,fcri-2; Wn-l,t^«-2)----Alfc^^(fc2, ^1 + i^; ^2, ^1 + ^) 

(^)From a more general point of view, one can obtain all the MCT predictions using a Landau expansion 
of F(G) in G{t,T) — G{t = oo, Tc), justified in dimension larger than six. This derivation, which unveils the 
generality of the MCT predictions, will the subject of a future publication. 
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where the matrix M is defined as: 

KA r ^ rr2R{ki,'^i)R{k-ki,ijj-ijJi) - - - - 

■^fe,a>('^2,fci;a;2,wi) = p{k- kiY — ^^^'^ ~ k\,u}2 - uii)V{ki,k2,k), 

where p is the particle density and C is the correlation function and R the response function 
(assumed to be related to C by the fluctuation-dissipation theorem (^)) for the density in 
Fourier space, and 

V = {ki- [hcih) + {h - k2)c{ki - fc2)]}{(fci - k) ■ [02 - k)c{k2 -k) + (fci - k2)c{ki - h)]} 

where c(fc) is the direct correlation function. If one follows the so-called schematic approxima- 
tion [11,12] where all structure factors S{k) appearing inside integrals are approximated by 
AS{k — ko) and all ingoing and outgoing momenta have modulus ko, then all fc-dependence dis- 
appears. Replacing the resulting (fc-indcpcndcnt) four leg vertex T'^V = T^S{ko)koA'^/{ATr'^p) 
with T^Vs-spin = 3 (where T is the temperature), one gets exactly the same diagrams of the 
mean field p = 3 case. For T <Tc and for small frequencies, the correlation function acquires 
a (non-crgodic) contribution: C{u)) = fS{uj) + Crcgi^) (where / is usually denoted q in the 
spin-glass literature) , while the response function tends to a constant at low frequencies given 
by (1 — f )/T (with power-law corrections in w, see below). Therefore, one obtains, as the dom- 
inant contribution: A„X4(w,f^) = C{^)5{oj)f\V{l - ff f]^'^-^\ Clearly, the series diverges 
when (1 — fYfV = 1, an equation precisely satisfied at the mode coupling temperature Tc, 
and signaling criticality (see, e.g. [11,12,15]). Hence, the asymptotic value of X4(t oo) is 
divergent at Tc- One can show, using a transfer matrix method to analyze all diagrams gener- 
ated from the eight building blocks drawn in Fig. 1, that the above singularity is unchanged, 
i.e. X4 indeed diverges as (1 — (1 — f^fV)^^ (^). Since / — jc^ \ft for T < Tc, we conclude 
that X4 diverges as e~-^/^ for T / T~ , indeed in agreement with the results of [14]. 

Our aim in the following is to understand in details how this divergence is affected by non 
zero frequencies and wave-vectors. In order to do so, we rely on the detailed results known 
about C and R in the vicinity of T^. One knows for example that when T = T,, — e, the plateau 
in C is reached after a time rj ^ e"^/^", where a is a non-trivial exponent. In the regime 
ry^ <^ u; <C 1, the response function acquires an extra contribution proportional to [11]. 
On the other hand, when T = Tc + e, two time scales diverge with different exponents. One 
is again the time r/ to reach the (pseudo-)plateau value / of the correlation function, and the 
second is the terminal time rj ^ e~'^ (with 7 = l/2a -I- 1/2&), beyond which the correlation 
finally drops to zero. Both the exponents a and b are non universal; for example, in the 
mean-field p = 3 case one has a w 0.395 and 6=1. Using these results for T < Tc and for 
non zero ui (corresponding to the shift in the measurement times for the four-point correlation 
function) we find that the term 1 — (1 — f)'^fV ~ ^/e leading to the singularity, is replaced 
by + Zlo'^, where Z is a constant. Thus, the only characteristic time scale is, for T < Tc, 
the plateau time r/. Resumming the series A„X4 is much more subtle in the case T > Tc. 
A naive transposition of the results for T < Tc to T > Tc suggests the following: since the 
l/y/e divergence of X4 for T < Tc can be traced back to /c — / ~ y/e, then the divergence 
for T > Tc should be 1/e because now f — fc — e. Furthermore, from the form of the ladder 
diagrams one expects that both the correlation time scale in the t direction and the peak time 

(^)This property is certainly true for the p = 3 spin-glass for T > Tc but may be problematic for MOT. 
However, we believe that the following conclusions are independent of the strict validity of fdt. 

(^)Actually a complete calculation should also take into ax;count the diagrams with extra vertices coming from 
imposing an initial condition at equilibrium. Although these diagrams do not change the singular behavior 
for T — » oo they are responsible for the vanishing of X4 when t oo. 
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Fig. 2 - Four-point susceptibility Xii''') in mean-field p = 3 spin-glass, obtained by integrating the 
dynamical equations obtained in [14], for various temperatures T > Tc {Tc = 0.612). The inset shows, 
in log-log scale, the position t* (black circles) and the height X4(**) (black squares) of the peak as a 
function of T — Tc, and the plain lines are power-laws with exponents —1 and —(1 -|- a) /2a = —1.76, 
corresponding to the 'naive' argument discussed in the text. We confirmed these results by also 
solving numerically the mixed p = 3 -f 6 case, for which a = 0.345 and h — 0.717 (open symbols). 
These last curves have been shifted for clarity. 

t* of Xii''') given by the terminal time scale tj. The above predictions turn out to be in 
perfect agreement with the numerical solution of the p = 3 mean-field case shovifn in Fig. 2 
(^), and can in fact be obtained analytically using another route [19]. 

Returning to mct in finite dimensions, one has to include the wavevector dependence of 
C, R and V. Assuming ergodicity is broken, Xi for fc = is found to diverge when the largest 
eigenvalue of matrix M is equal to unity: 

M(fc2,fci) = ^S{h)S{k2)S{ki - k2)V{kuk2,0Kl - (4) 

with, for T < T,: /,/(! - /,) = J d^k' /{2nfS{q)S{k')S{q~ k')V{q,k' ,Q)hf^_p. The 
spectrum of M has been previously studied in detail because it controls the convergence of 
the above equation on fq when solved by an iterative procedure. It has been established 
(see [12]) that the maximum eigenvalue of M is non degenerate and approaches one as ^/e at 
the transition, which in turn leads to the famous mct singularity of f — fc- The effect of a 
non zero wave- vector k can be understood using perturbation theory. By symmetry, it is clear 
that the correction to the largest eigenvalue must be, for small fc, of order k^. This leads for 
small k and T < Tc to a propagator (Ffc^ -I- y/e)^^. Naively transposing the above results for 
T > Tc gives a propagator behaving as (F'fc^ -|- e)~^. Thus, we conclude that the four-point 
correlation computed, say, for t = and S of the order of the particle radius, and t of the order 
of the correlation time scale (r/ for T < Tc and Tf for T > Tc) will behave as the two-point 

(*)Note that these numerical results do not coincide with the results reported in [14], obtained with the same 
code; in particular, the height of the peak was found to diverge as e~^/^ and not as e~^. A possibility is that 
the 'conjugated field' used in [14] was not small enough. 
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correlation in standard critical phenomena, i.e as G (j) /r'^~'^+^. From the above results, we 
find ?7 = 0, and a dynamical length (, that diverges as e^" , with v = 1/4, 1/2 respectively 
below and above the transition. The above scaling law for T > is compatible with recent 
numerical simulations [5-7]. It is also interesting to note that for T = Tc and for small k,co 
our results imply that ^ behaves as Ffc^ + Zlu°', which allows us to identify the dynamical 
exponent as z = 2/a. At the critical point, length and time scales are related by an anomalous 
sub-diffusion exponent: r ^ i"/^, with a < 1/2 [12]. The above result will hold for T > Tc 
and tJ^ <C w <C 1, whereas for UTt <C 1 the denominator of the propagator will rather behave 
as + (Ttui)'^, corresponding to simple diffusion on long time scales. Note that the relation 
between £ and allows one to define a second dynamical exponent z' = l/a+1/6 — 2j, echoing 
the presence of two diverging time scales for T > Tc. Numerically, for the Kob- Andersen LJ 
mixture, z' = 4.5±0.2 [7], a value surprisingly close to our prediction 27 = 4.68 [12]. However, 
exactly as for usual critical phenomena, one should expect long-wavelength fluctuations to be 
dominant below some upper critical dimension dc and changing the value of all the exponents, 
at least sujjiciently close to Tc- The value of dc can be obtained by analyzing the diagrams 
correcting the mct contribution to the self-energy considered here, or using a Ginzburg-like 
criterion; both lead to dc = 6. The diagrammatic is in fact very similar to that of the (j>^ 
theory, which is in a sense expected since the order parameter is here the correlation function 
itself, which does not have the Ising symmetry. More physically, one can argue that for the 
spatial fluctuations to bo irrelevant, these should not blur the \/e singularity for T < Tc (or 
the e singularity for T > Tc) of the non-ergodic parameter / around fc. Within a sphere of 

radius these fluctuations are of order i"^ J r'^^'^Gi{r)d.r ^ to be compared to 

the total contribution of the singularity, I'^^e for T <Tc (or I'^e for T > Tc). Using e ~ i''^ 
for T < Tc (e ~ £~'^ for T > Tc), we thus find that the fluctuations become dominant for 
large £ whenever {d + 2) /2 > d — 2, or d < dc = 6. A detailed calculation of the exponents for 
d <6 would be very interesting, to estimate how the mean field exponents of MCT are affected 
by spatial fiuctuations in d = 3. A naive guess, based on a percolation interpretation of the 
MCT transition [20], suggests z/ ~ 0.88 and i] ~ 0. [The analogy with phase transitions could 
also shed light on the non trivial (fractal) structure of the mobile regions, see [21-23].] The 
knowledge of these critical exponents is crucial, since a quantitative fit of many experimental 
results has been attempted using MCT [12]. In particular, neither the ^/f_ singularity of the 
non-ergodic parameter for T < Tc (argued to be a signature of the MCT singularity) nor the 
exponents derived above for X4 are expected to remain valid for = 3. Nevertheless the 
existence of a diverging length scale could hint at a large (and perhaps unexpected) degree of 
universality in the dynamics of glassy systems [5,7]. 

The critical fiuctuations discussed above should however not be confused with another type 
of fluctuations that are expected to destroy the MCT singiflarity altogether, and suppress 
the glass transition. These fluctuations are the 'activated events' or 'hopping processes' that 
prevent a complete freezing of the super-cooled liquid, but the consistent inclusion of these 
in an extended version of MCT is still quite a challenge [12, 16]. The way the two types of 
fluctuations interact to make any of the above flnite dimensional MCT predictions observable 
is unclear to us: a second Ginzburg criterion, pertaining to activated fluctuations, should be 
devised and compared to tlw^ one abov(\ 

The existence of a diverging dynamical correlation length puts strong constraints on any 
theory of glass forming liquids. We have shown that finite dimensional MCT does, at least 
qualitatively, survive the test (note that the existence of a diverging length scale in MCT 
should allow for the observed decoupling between diffusion and viscosity.) Only a more quan- 
titative comparison of the predictions of MCT in three dimensions would allow one to rule it 
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out entirely, or to confirm that it is indeed a useful picture (at least sufBciently far from Tc). 
Simulations of systems with different fragilities and/or in higher dimensions would also be 
very interesting to gain further insight, and test different predictions or scenarii. For instance, 
the critical mobility defect scenario of [7] predicts an upper critical dimension dc = 4 with the 
same dynamical exponent z w 3.7 for both strong and fragile-to-strong liquids, in contrast 
with the prediction of spatial mct, where dc = 6 and where both a and b (and therefore z) are 
system-dependent. Finally, the extension of the ideas proposed here to non-zero shear rates 
and aging situations would be extremely valuable (see [23,24]). 

* * * 
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